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Abstract 

|Zj ■ The ground state of the ring-shape magnetic nanoparticle is studied. Depending 

on the geometrical and magnetic parameters of the nanoring, there exist differ- 
ent magnetisation configurations (magnetic phases): two phases with homogeneous 
magnetisation (easy-axis and easy-plane phases) and two inhomogeneous (planar 
^ . vortex phase and out-of-plane one). The existence of a new intermediate out-of- 

OO I plane vortex phase, where the inner magnetisation is not strongly parallel to the 

*vg ' easy axis, is predicted. Possible transitions between different phases are analysed 

Tij" ■ using the combination of analytical calculations and micromagnetic simulations. 
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1 Introduction 



Magnetic nanoparticles and their structures have become a subject of interest 
in the past few years [1[. For nanoparticles of axial geometry, such as disks 
(nanodots) and rings (nanorings) , the vortex configuration becomes the low- 
est energy state when the particle diameter exceeds the single-domain size 
owing to the competition between exchange and magnetostatic interactions 
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J2]. Because of nontrivial topological properties, the vortex state nanoparti- 
cles are promising candidates for the high-density magnetic storage devices, 
some designs of magnetic memories and high-resolution magnetic field sensors 



There are known two types of vortices in a nanomagnetism. For the disk- 
shaped particle there exist an out-of-plane vortex (OPV), which is charac- 
terised by the core magnetisation perpendicular to the disk plane. Different 
static and dynamical properties of magnetic nanodisks have been studied re- 
cently. In particular, theoretical and experimental studies of the transition 
between single domain and vortex configurations (phases) [3|, y, |5|, la, LD] give a 
possibility to conclude, that the vortex state in a nanodisk can exists for the 
dot diameter in the submicron range, which is much more than the minimum 
nanoparticle size, currently used in experiments. That makes serious limita- 
tions for considering such structure as competitive one for the magnetisation 
storage. 

The ring shape potentiates more magnetisation configurations, see for review 
J8|. It is well-known [8j that the pure in-plane vortex (IPV) configuration can 
be realised in the nanoring, where the perpendicular magnetisation component 
is absent. Such a vortex has lower energy than the OPV in a disc, it does 
not produce surface magnetostatic charges, thus it can provides more stable 
magnetic bit for the magnetisation storage. 

Since there exist pure OPV for the disk and pure IPV for the ring, there 
appears a reasonable question, whether the intermediate state exists or not. We 
will show in the paper that intermediate vortices (IMV) exist in the nanoring 
when the inner hole is small enough. For such vortices the magnetisation inside 
the vortex core has an out-of-plane component, but even in the vortex centre 
it is not perpendicular to the ring plane. The presence of the hole in the disk 
decreases the "singularity" in the vortex centre and so it gives an opportunity 
to realise the vortex state in a smaller nanoparticle. We will demonstrate that 
the vortex state can be realised in a nanoring with diameter of few tens of 
nanometers. 

In this paper we provide the systematic study of the magnetic nanoring's 
ground state. Namely, by combining analytical methods and micromagnetic 
simulations we describe the equilibrium magnetic phases: two homogeneous 
states with magnetisation directions along the ring plane and perpendicular 
to it, and inhomogeneous vortex states, including OPV, IPV, and IMV states. 
The key moment is to propose the simple analytical approach, which allows 
us to describe different states and transitions between them. 

The paper is organised as follows. In Sec. 2 we formulate the model and treat 
it analytically for the homogeneously magnetised nanoring. The vortex-state 



nanorings are analysed in Sec. 3 both analytically, using simple two-parameter 
ansatz, and numerically, using micromagnetic simulations. The phase diagram 
of possible ground states compares analytical predictions and numerical sim- 
ulations (Sec. 4). 



2 The model. Homogeneous states 



We start with the continuum description of the spin distributions inside the 
nanoring in terms of the magnetisation, normalised by its saturation value, 
m = M/Mg, and content ourselves with two main contributions to the energy 
functional, the exchange energy Ecx and the magnetostatic energy Ems- We 
neglect the anisotropy energy contributions, what is a reasonable for the soft 
materials like permalloy (Ni8oFe2o, Py)- 

For the homogeneous magnetisation distribution the total energy contains only 
the contribution of surface magnetostatic charges: 

Iff {M . dS){M' . dS') 



Let us consider the nanoring, R is the outer radius, a is its inner radius, and 
h is the thickness. First of all we study the case of the ring homogeneously 
magnetised along its plane. It is convenient to consider the energy density, 
normalised by the value irR'^hM^: 

^ , , pdxdCp'dx'dC cos X cos x' 

vVms 
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where S indicates the outer and inner edge surfaces, (p, x, C) ^'^^ the cylindrical 
coordinates, p is the radius normalised by the magnitude of the outer ring 
radius i?, the variable (^ is the thickness normalised by h. We have introduced 
also the particle aspect ratio e = h/2R, and radii ratio a = a/R. Integration 
over edge surface can be rewritten in explicit form: 

W^s = '^(l. 1) - Sa^^ll, a) + tt^(?(«, «), (2a) 
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Straightforward but long and tedious calculations result 
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Here K(a;), E(x) and Il{x,y) are the complete elliptic integrals of the first, 
second and third kind respectively. 

In the limit case of the disk, the bulk expression (2) transfers into the well- 
known m formula 
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To calculate the energy of the ring, homogeneously magnetised along z-axis, 
one can use relations between the magnetometric demagnetisation factors (lfll |. 
These may be written in the form: 



2W^MS + Wms 
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(4) 



where the factor (1 — a^) appears in consequence of our normalisation of the 
energy density. 

It is well known that the homogeneously in-plane magnetised (easy-plane, EP) 
state is preferable energetically for thin enough disk particles; when the disk 
aspect ratio is greater then the critical value Eq ~ 0.906 [9], the homogeneously 
out-of-plane magnetised (easy-axis, EA) state is realised. 

Here we make a similar critical analysis for the nanoring. At critical values 
of the parameters W^^ = W{^^, and the critical curve eda) can be found by 
solving the equation W^g{ec, a) = ^(1 — a^). The numerical solution of this 
equation is plotted in Fig. 1. For the approximate description one can use the 
asymptotically correct solution 
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Fig. 1. The critical curve for the homogeneously magnetised states. Symbols corre- 
sponds to the numerical solution (see the text), and the curve is the approximate 
solution (5). 

which reproduces the numerical results with the accuracy within 6 x 10~^, see 
Fig. 1. 



3 Vortex state 



The phase diagram, see Fig. 1, describes the single-domain formation only. 



When the particle size exceeds the exchange length, Zgx = v ^/47rM^ (A is the 
exchange constant of the material), a magnetisation curling takes place [2], and 
the vortex-state can be energetically preferable. We start with the exchange 
energy in the form E'ex = 2^h J d^a:( Vm)^. Using the angular parametrisation 
for the normalised magnetisation m = {cos ^ cos 0; cos ^ sin 0; cos ^}, one can 
describe the vortex solution as follows: 



9 = 9{r), 
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Here (r, x) are the polar coordinates in the ring plane. For the vortex-like 
solution (6) the exchange energy density takes a form: 
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Let us calculate the magnetostatic energy. For the vortex distribution (6) 
the volume magnetostatic charges are absent (V ■ m = 0), and the surface 



contribution (1) can be presented in the form J2| 
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At this point we should specify the profile function 9{p). There are known sev- 
eral models for describing the vortex structure in the disk-shape particles. We 
should mention the core model by Usov and Peschany [ll| and its modifica- 
tions |a,ll2,ll35 which neglects the magnetostatic interaction inside the vortex 
core, and also has a singularity at the core radius. More realistic model was 
proposed by HoUinger et al. p , where the vortex has a bell-shaped structure. 
This model is in a good agreement with simulations data for a disk. Neverthe- 
less we will use another model by the following reasons: (i) the model should 
be able to describe magnetisation distributions both in disks and rings, and 
(ii) it should be very simple to be treated analytically. All mentioned above 
models can not be simply generalised for the ring case, and can not be analysed 
analytically. 

In pursuance of above mentioned reasons we propose the following two-parameters 

ansatz: 

r 
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Here parameter /i describes the vortex amplitude in the centre of the ring. 
The typical vortex width is determined by the parameter A. An advantage of 
this ansatz is a possibility to describe different kinds of vortices: the OPV for 
fi = 1, the IPV for /i = 0, and the IMV for fi G (0, 1). For the case of the disk 
(yU = 1) our ansatz (9) fits well the structure of the pure OPV in easy-plane 
magnets jl4| . 

For the analytical treatment of the model, we make also one serious simplifi- 
cation: we use the local shape-anisotropy model [15] instead of the nonlocal 
magnetostatic energy (8), what is acceptably for thin particles (e: ^ 1): 

H^anis = -^y^ rdrcos2^(r). (10) 



Finally the total energy density for the vortex state ring reads 
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This integral can be derived analytically, the explicit form is calculated in Ap- 
pendix A. We have derived the energy (11) using the local shape-anisotropical 
model (10) instead of the non-local magnetostatic one (8). Strictly speaking. 
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Fig. 2. The comparison of two OPV profiles : the ansatz (9) with // = 1 and 
A = v2 (sohd red hne) and micromagnetic simulations (black bars) for the Py 
disk (i? = 50nm, h = 5nm). 

the simplified model is valid only for infinitesimally thin rings, see Appendix 
B for details. However further analysis demonstrates that it provides a rea- 
sonably accurate picture of phenomenon even for a large aspect ratio. 



3.1 Pure out-of-plane vortex 



We start our analysis with the case of a disk-shape nanoparticle, where the 
pure OPV is realised. This well-known case is a good test for our simple 
theory. The magnetisation in the vortex centre is perpendicular to the disk 
plane, therefore the vortex amplitude at origin is equal to unit, /i = 1. 

The energy of the vortex state nanodisk can be easily derived as a limit case 
of Eq. (11), see Appendix A for details: 
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Here ^r = 2R'^ / X^l"^^, 7 ~ 0.577 is Euler's constant; Ei(^) is the exponential 
integral function. The value of the variational parameter A can be found by the 
minimisation of the energy (12) with respect to A. Finally, the A-parameter 
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Fig. 3. Micromagnetic simulation data of the IMV profiles for different radii of inner 
hole. Other parameters are the same as in Fig. 2. 

can be calculated as a solution of transcendental equation 
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which can be solved numerically. For the case of large enough disk radii {R > 
4/ex), the numerical solution A(-R) practically coincides (with the accuracy 
within 2 x 10^^) with the limit value A(oo) = y^- 

To verify our analytical model calculations we have preformed the numeri- 
cal computer simulations using a three-dimensional OOMMF micromagnetic 
simulator code 16|. In all simulations we have used the following material pa- 
rameters for the Py: A = 1.3 x 10~® erg/cm (using SI units A^^ = 1.3 x 10^^^ 
J/m), Ms = 8.6 X 10^ G (Mf^ = 8.6 x 10^ A/m), and the anisotropy have been 

neglected. This corresponds to the exchange length /ex = yA/AirMg ^ 5.3nm 



'A/hqMs). a comparison of the vortex profiles, obtained by simula- 
tions and by analytical approach is presented in Fig. 2. Our ansatz (9) fits the 
simulation data with the accuracy within 0.09. 



3.2 Intermediate vortex solution 



Let us consider the ring-shape nanoparticle. The presence of the hole in the 
centre of the disk changes the topological properties of the vortex solution. 
Hence it is not necessary for the vortex amplitude to be equal to unit in the 
ring centre. Moreover, if the inner hole is big enough, there can exist pure 
IPV, where the magnetisation does not have out-of-plane component at all. 
In the intermediate case of small enough inner hole the vortex amplitude /i 
should vary in the range [0, 1]. We have verified this idea by the micromagnetic 
simulations, as described above. One can see from Fig. 3 that the out-of-plane 
vortex structure has a well-defined bell shape, but the amphtude of the vortex 
decays fast when the inner hole becomes bigger. When the inner hole is greater 
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Fig. 4. Vortex amplitude (at the inner edge) vs inner radius of the ring. Bars denotes 
simulations results and solid green line corresponds to theoretical calculation, which 
were performed numerically. Dashed blue line corresponds to the fit (16). Other 
parameters are the same as in Fig. 2. 



than some critical value Oc, which is about 2nni for the Py ring, only pure IPV 
solutions can be realised. 

We have chosen the ansatz (9) because the bell shapes of the vortex structures 
are well pronounced for different vortex amplitudes. Using this ansatz we 
have derived the general expression for the vortex energy (11). To calculate 
variational parameters of our model /z and A, we have minimised the energy 
(9) with respect to these parameters: dW/dfi = and dW/dX = 0. This set 
of equations can be solved only numerically. Analysis shows that the solution 
practically (with the accuracy within 1 x 10~^) does not depend on R when 
R > 3/ex, hence one can use the limit case i? ^ cxd in calculations. The results 
are presented in Fig. 4. Both analytical and simulation results have the same 
main behaviour: the vortex amplitude slightly decays when the inner hole 
increases for small enough holes, but it sharply breaks at some critical inner 
radius Oc- However, the value of the critical inner radius is about 2 nm from 
simulations while Oc = 1-6 nm from the theory (see below). This discrepancy in 
critical values is caused by the local model of the magnetostatic interaction. In 
Appendix B we will take into account the nonlocal magnetostatic interaction, 
which gives the critical inner radius 1.94 nm, what is in a good agreement 
with simulations results. 



To calculate the critical value ac, when the transition to IPV occurs, one can 
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Fig. 5. The numerical solution of the Eq. (15). 
expand the IMV energy density (11) in series on /i at /i = 0: 
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One can see that C2 is always positive, hence the vortex solution with /i 7^ 
can exist only when Ci < 0. This is the case of the double-well potential with 
minima at /i^ = \Ci\/C2- The pure IPV solution with /i = corresponds to 
the case of Ci > 0. Thus the critical point can be found from the condition 
Ci = 0. The numerical solution of the equation 
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is shown in Fig. 5. It is not necessary to check the second minimum condi- 
tion with respect to A. Obviously, the critical parameters correspond to the 
maximum of the curve a(A). Finally, the critical inner radius a^ ~ O.S/cx, see 
Appendix A. 

At vicinity of the critical point, the vortex amplitude has asymptotic /x oc 
y^oT^^^- One can extend this result and use the function 



/i 



approx 



1 — a/ttc 



(16) 



which fits the numerical results of Fig. 4 with high accuracy (about 9 x 10 ^) 
in the whole range of the parameters. 

When the inner ring radius a exceeds the critical value Qc-, the pure IPV 
solution with /x = takes place. In this case all magnetostatic charges are 
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Fig. 6. Three-dimensional phase diagram of magnetisation ground states from an- 
alytical analysis. 

absent, and the vortex energy density has only exchange contribution, which 
simply reads 



W^ 



IPV 



4<, R 



(17) 



4 Phase diagrams 



Summarising results on the equilibrium magnetisation distribution, we have 
calculated energetically preferable states for different ring geometries. Since we 
have three parameters, which define the ring, i?, a, an /i, our phase diagram is 
a three-dimensional plot. Theoretically calculated phase diagram is presented 
in Fig. 6. Different phases are separated by boundary surfaces. The general 
properties of the phase diagram are as follows. The ground state of the very 
thin ring is the homogeneous easy-plane state. At the thickness increasing we 
can switch either to the homogeneous easy-axis state, when the outer ring 
radius is small enough, or to the vortex state, when the ring is large. The 
boundary surface between two homogeneous states can be simply expressed 
analytically, h = 2Rec{a), where the £:c(a)^dependence is described by Eq. (5), 
see also Fig. 1. Now, if we increase the inner radius of the ring, the vortex 
state becomes energetically preferable. One can follow step-by-step transitions 
between different vortex phases: from the pure OPV state for the case of the 
disk (a = 0) to the IMV state for the ring with small radii ratio (a < Oc) and, 
finally, to the pure IPV state when a > ac- 

To verify our theoretically calculated phase diagram, we use 3D micromagnetic 
simulations, as described above. Simulations data are reproduced in Fig. 7 to- 
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Fig. 7. The phase diagram of magnetic ground states for different radii ratio. The 
areas of different colours correspond to different kinds of magnetisation distribution. 
Symbols represent the simulation data: squares - EP, diamonds - EA, disks - IMV, 
disks with vertical lines ~ OPV, disks with horizontal lines - IPV. The dashed 
lines in fig. (a) and (b) are the theoretical borders calculated within the nonlocal 
model of magnetostatic interaction. Solid lines correspond to the theoretical borders, 
calculated within the local model. 



gether with theoretical results. Fig. 7a presents the phase diagram for the case 
of the nanodisk. Numerical results for this case is in a good agreement with 
previous theoretical ja |l2| and experimental J7| results. One can compare the 
theoretical results obtained within the local (solid line) and nonlocal (dashed 
line) models of the magnetostatic interaction. In the case of thin disks the 
boundaries are coincident and agree with simulations results. However in the 
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Fig. 8. The triple point analysis: the outer ring diameter as a function of (a) the ring 
thickness and (b) the radii ratio. Red lines corresponds to the exact numerically 
calculated data (the local model of the magnetostatics) , dashed blue line is the 
analytical solution (18), which is acceptable for a > ao. Points with error bars were 
obtained from simulations data (Fig. 7) 

case of thick disks for the boundary between OPV phase and homogeneous EA 
solid line (the local model) lies higher then simulations data. This is because 
the variational approach in the local shape-anisotropy model always gives the 
upper limit of the energy. 

The presence of the hole inside the disk drastically changes the vortex state. 
The vortex amplitude depends on the inner radius of the ring. If the radii ra- 
tio is small enough, the phase diagram contains four phases, see Fig. 7b. The 
boundary between phases IMV and IPV was calculated using nonlocal model 
of magnetostatic interaction (see Appendix B). We have not managed to cal- 
culate the boundaries between homogenous and vortex states within nonlocal 
model, because in the case of ring there are two variational parameters (/i and 
A), in contrast to the case of disk, where only one variational parameter (A) is 
present. If the inner radius exceeds the critical value Qc, only one vortex phase 
with IPVs can exist, see Fig. 7c and Fig. 7d. 

The triple point in the phase diagrams is the object of a special interest. In 
particular, it provides an information about the lowest possible radius of the 
vortex state nanoparticle. The triple point was analysed for the disk particles 
using the scaling technique 13|,|l7|. Here we present results for the case of the 
ring. Since at this point the vortex energy Wiuy is equal to the homogeneous 
state energies, moreover W^^ = W^^, one can calculate the triple point pa- 
rameters by solving an equation VFimv = ^(1 ~ o;^)- ^^ general case it can 
be done numerically only, see Fig. 8. An analytical analysis can be performed 
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asymptotically for small a, where we can limit ourselves to the linear depen- 
dence 2R^^^'> = h^^^^ /eq, where Eq ~ 0.906 is the critical aspect ratio for a disk 
J9|. In the limit case of a = the triple point is characterised by the following 
parameters: i?^*''-' ~ 3.73Zex and h^^^^ ~ 6.76/ex (The values were obtained for 
the nonlocal model of the magnetostatics). By increasing the inner radius, the 
triple point radius R^^^^ decreases, so the IMV can exist in smaller particles 
than the pure OPV. At some critical inner radius Oc the transition to the IPV 
state occurs. The energy of the pure IPV has a simple form (17), so the triple 
point analysis can be easily done for this case: 
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a^ a 



where ec(a)^dependence is described by Eq. (5). Such simple analytical de- 
scription of the triple point is valid only, when a > ac and the IPV is stable. 
Simple calculations show that it can be realised for a > ao ~ 0.077. The triple 
point radius i?^*''-' and thickness /i^*''-' diminish when a increases, and they have 
the asymptotic behaviour: 

R^''\a) ~ l^^V3\l + {l-a)/2\, h^''\a) ~ 2k,V3{l-a). (19) 



a^l 



a^l 



The limit value R^^^\l) = IcxV^ provides the lowest bound for the vortex state 
magnetic ring. For the Py nanoring this limit is about 9.18nm. 



5 Conclusions 



We have presented a detailed study of the ground state of magnetic nanorings, 
including homogeneous (easy-axis and easy-plane) states and inhomogeneous 
vortex states. In addition to conventional out-of-plane and in-plane vortex 
types we have found a new type of vortices in a nanoring. Such intermediate 
vortices are characterised by the smaller amplitude in the vortex centre, they 
are the key point to understand the transition from OPV to IPV state. We 
have proposed also a simple analytical description of the IMV state. Using 
this approach we have studied analytically phase transitions between differ- 
ent ground states. All results are confirmed by our direct 3D micromagnetic 
simulations. 

Using the critical analysis of the phase diagram triple point, we conclude that 
the lower bound of the vortex-state nanoparticle decreases for a ring geometry 
in comparison with a disk. The minimal nanoring radius varies from about 
3.7/ex for the case of the disk to V^hx for the case of narrow ring. 

Our model of the IMV can be applied for the description of dynamical effects 
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of the switching of the vortex polarisation [18|, |l9|, |20|. We expect that the 
dynamical switching effects should be more pronounced in the ring geometry, 
because the vortex amplitude /i can smoothly vary and even change its sign. 
This is the subject of future research. 
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Appendix A The vortex state energy calculations 



Let us start with the vortex energy in the form (11). Using the direct integra- 
tion one can write down the energy in the form 
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Here ^a = ,r, and C,r = ; the integral I{^) can be expressed as follows 
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where Ei(a;) is the exponential integral function. 
Using the asymptotical behaviour 

m ^~, ^ + /^^ - In e(l - /i^) + / ^^^^dt, (A.3) 
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one can easy calculate that X(0) = A^/2 + 7 when /i = 1. According to (A.l) 
the vortex energy of a disk with radius R can be written down as 



W< 



OPV 



2vr^ 
R^ 



no) - l{in) 



^Ji=l 



(A.4) 



which results in (12). 



Let us go back to the case of the ring and calculate the critical inner radius, 
when the transition from the intermediate vortex solution to the pure in-plane 
one takes place. One can rewrite Eq. (15) as follows: 



F(a,A) 



dx xe~ 






0, 



:i5') 



which is an implicit form of the dependence a(A). The critical value Oc can be 



calculated from the condition dF/d\ = 0; hence Oc = IcxJxo — Xq, where xq 
is a positive root of the equation 2e~^ + Ei(— x) = 0. This results in xq ~ 0.1 



and a^ ~ 0.3/cx- 



Appendix B Magnetostatic energy calculations 



Let us write down an expression for the vortex energy with account of the 
non-local magnetostatic contribution (8): 



W, 



IMV 



-^ I ^dx 






+ A^/i^/C X, 



/ 


4x2 


V- 


- /i^e"^^^ 


a 


R h 



x^ 



^'ex ^^ex ^'ex , 



(B.l) 



Here the kernel of integral 

-1-2 

]C(x,C,ri,0 = —- x'e-'' / (l-e-'^*)Jo(xt)Jo(x't)dtdx'. (B.2) 
C Ji JO 

Note that the kernel of the local model (11) has an exponential decay, /C(x, •, •, 
exp(— 2x2). 

Let us analyse the case of small thickness {( ^ 0): 



r.12 



/C(x,e,^,C ^ 0) = e"^' r x'e-""'^ I tJo(xt) Jo(x't) dt dx' 

Ji Jo 



(B.3) 
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The internal integral in (B.3) can be calculated directly: 

oo 

f tJo{at)Jo{bt)dt = lim — — - [aJo{bt)Ji{at) - 6Jo(at) Ji(6t)] = 

(B.4) 

1 j sin(a - b)t cos{a + b)t \ 1 , n ^ -.. ^ ^i 

*^°° Vab [ 7r(a - 6) n{a + o) J Va6 

Finally the kernel (B.3) can be written in the form 

/C (x, e, r], 0) = e-"' r J- e"""5(a:' - x)dx' = e^'"'. (B.5) 

Thus the nonlocal model of magnetostatic interaction passes into local model 
in the limit of thin films. 

An account of nonlocal magnetostatic interaction changes the value of the 
critical inner radius, when the magnetic phase transition into the pure IPV 
phase occurs. This critical radius can be calculated numerically as a solution 
of the following equation: 



oo 

dx X 



X / \ Alp^ \i, 



e-^-UAx^-^\+\^lC\x,^,oo, 



ex , 



(B.6) 



To simplify the numerical calculations one can reduce the kernel to 

;C(x, e, oo, C) =^| r ^^:!^ii^dt + r e-*^K {x/t) dt 

fB 71 

A v'^^ + (^ + ty 

where K(x) is the elliptical integral of the first kind. 

The solution of the Eq. (B.6) with kernel (B.7) is shown in Fig. B.l. The 
maximum of dependence a{\) that corresponds to the critical value of inner 
radius is Oc ~ 0.37/ex- For the Py nanoring ac ~ 1.94 nm, which agrees with 
the simulations data, see Fig. 4. 

An account of nonlocal interaction gives a possibility to calculate the critical 
inner radius on thickness of the ring, see Fig. B.2. This dependence can be 
fitted by the function 



af(/i) ^aeV/l + cV/ex, c^O.68, (B.8) 
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Fig. B.l. The numerical solution of the Eq. B.6 is denoted by solid red line. During 
the calculation it was assumed that parameter — = ^^^^ , that corresponds to the 



parameters of the simulations (see Fig. 4). Dashed black line corresponds to the 
solution of the model of local magnetostatics (15), Fig. 5. 




5 10 15 

Ring thickness (/i//ex) 

Fig. B.2. The dependence of the critical inner radius on thickness of the ring. Red 
bars denotes the exact values, which were obtained numerically. Solid black line 
corresponds to the fit (see text). 

where Oc ~ 0.37/ex is the limit value for infinitesimally thin rings, which was 
calculated in Sec. 3. We also have numerically calculated the dependence of 
vortex width (A) on disk thickness (for nonlocal magnetostatic model), see 
Fig. B.3. 
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Fig. B.3. The dependence of the vortex width in case of disk on disk thickness. 



This dependence also can be well fitted by the third-root function: 



\^\h)^V2m+dh/h^, t/^0.39. (B.9) 
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